arXiv: 1501.06345v2 [astro-ph.CO] 18 May 2015 


Oscillating modulation to B-mode polarization from varying 
propagating speed of primordial gravitational waves 

Yong Cai 1 ’*, Yu-Tong Wang 1 '*', and Yun-Song Piao 1,2 * 

1 School of Physics, University of Chinese Academy of Sciences, Beijing 100049, China and 
2 State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, P.O. Box 2735, Beijing 100190, China 

Abstract 

In low-energy effective string theory and modified gravity theories, the propagating speed ct 
of primordial gravitational waves may deviate from unity. We find that the step-like variation of 
ct during slow-roll inflation may result in an oscillating modulation to the B-mode polarization 
spectrum, which can hardly be imitated by adjusting other cosmological parameters, and the 
intensity of the modulation is determined by the dynamics of ct ■ Thus provided that the foreground 
contribution is under control, high-precision CMB polarization observations will be able to put tight 
constraint on the variation of ct , and so the corresponding theories. 
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I. INTRODUCTION 


Inflation, as the paradigm of the early universe, has not only solved a lot of fine-tuning 
problems of the big bang theory, but also predicted the primordial scalar and tensor per¬ 
turbations. The primordial tensor perturbations, i.e. primordial gravitational waves (GWs) 
[1] [2] [3] , have arisen great attentions after the BICEP2 collaboration’s announcement of the 
detection of B-rnode signal in the CMB (around l ~ 80) [4], which was interpreted by them 
as the imprint of the primordial GWs, though this result is doubtful due to the foregrounds 
of polarized dust emissions [5] [6] , see also [7] . 

The detection of primordial GWs would verify general relativity (GR) and strengthen our 
confidence in inflation and quantum gravity [8], and also put more constrains on inflation 
models and modified gravity at the same time. Besides the CMB experiments which mainly 
aimed at detecting low frequencies (10 —1 ' — 10“ 15 Hz) GWs, many experiments relate to 
higher frequencies based on other methods, such as pulsar timing array (10~ 9 — 10 -8 Hz), 
laser interferometer detectors (1CT 4 — 10 4 Hz), will be carried out in the coming decades. 
However, since the amplitude of the GWs would stay constant after they are stretched 
outside the horizon, and decrease with the expansion of the universe after they reenter the 
horizon, the primordial GWs with longer wavelength provide the most of opportunities for 
the detection [9]. Therefore, the CMB observations, especially the CMB B-mode detections, 
are still the most promising experiments to detect the primordial GWs if they actually exist. 

Einstein’s GR is the most accepted theory of gravity. However, it might be required 
to modify when dealing with the inflation in the primordial universe and the accelerated 
expansion of the current universe. During the matter and radiation dominated era, modified 
gravity has several effects on the CMB spectra, such as the lensing contribution to B-modes 
[10] and the variation of propagating speed ct of primordial GWs [11][12], we are especially 
interested in the latter in this paper, see e.g.[13] for the case with the scalar perturbation. In 
GR, graviton is massless and propagates along the null geodesics, so the propagating speed 
ct of GWs is naturally set to be unity, i.e. the speed of light. But in modified gravity, e.g., 
the low-energy effective string theory with high-order corrections [14] [15] [16] [17] [18], and 
also modified Gauss-Bonnet gravity [19], and generalized Galileon (Horndeski theory [20]) 
[21][22][23][24], and beyond Horndeski theories [25],[26], and the effective theory of fluids 
at next-to-leading order in derivatives (e.g. [27]), ct might deviate from unity. Because the 


2 


value of ct determines the time of horizon crossing of GWs, during the recombination the 
change of cr can result in a shift of the peak position of the primordial B-modes, see [11] [12], 
for example. 

In this paper, we focus on the effect of the variation of Cr during slow-roll inflation on the 
CMB B-mode polarization, and show how it offers a distinct way to test the modified gravity 
theories. We find that the step-like variation of cr may result in an oscillating modulation 
to the B-mode polarization spectrum, which can hardly be imitated by adjusting other 
cosmological parameters. The intensity of the modulation is determined by the ratio of 
Cp before and after the variation, which depends on the dynamics of theoretical models. 
This oscillating modulation is so rich in feature that it may easily be discriminated from 
the variation of other parameters or other features. Thus provided that the foreground 
contribution is under control, high-precision CMB polarization observations will be able to 
put tight constraint on the variation of cp, and so the corresponding theories. 


II. OSCILLATING SPECTRUM OF PRIMORDIAL GWS 


We begin with the action for the GWs mode h i j 1 e.g. [IT] [23] 


S(2) 



l' 2 


,(V/i 


VJ 


( 1 ) 


where obeys d t h t .j = 0 and h tl = 0, Qt is regarded as effective Planck scale Mp e ^{r) 1 cp 
is the propagating speed of primordial GWs, and the prime is the derivative with respect 
to the conformal time r, dr = dt/a. During slow-roll inflation, the slow-roll parameter 
e = —H/H 2 <C 1, as well as 


e Q 


Qt «1, s =-?-«! 


( 2 ) 


HQ t Hc t 

are required, e.g., see [28] for a recent study. 

Here, we will mainly focus on the effect of varying cp, i.e., the condition s -C 1 might be 
broke at some point, on primordial GWs spectrum. Noting that the effects of varying sound 
speed of primordial scalar perturbations on the scalar power spectrum have been investigated 
in [29] [30]. We do not get entangled with the details of (1) and the evolution of background, 
and will assume that the background is the slow-roll inflation, which is not affected by the 
variation of cp, and set Qp constant and Mp = 1. We will discuss a possibility of such a 
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case in Appendix A. In addition, there may be a mass term [31] [32] in (1), which might also 
be time-dependent [33], but we will not involve it. 

We can expand into Fourier series as hj 3 (r, x) = f ^4-e _ * k ' x ^(T, k), where 


Kj (t, k) = ^ |M t ; A;)a A (k) + hi 


-k)a\ I 


k ) 




( 3 ) 


A—“I - , x 


where the polarization tensors (k) are dehned by kje^\ k) = 0, e^(k) = 0, which satisfy 
e-^(k)e*j A ' l (k) = S\y, e*j A ' ) (k) = e^(—k), the commutation relation for the annihilation 
and creation operators a A (k) and a^(k') is [a A (k), (k')] = <5 AA /5®(k — k'). We dehne 

h A (r, k ) = Uk{r )/zt and zt = ay/Qr/Z- Thus we get the equation of motion for u k as 

< + [elk 2 ~^)u k = 0. (4) 

To phenomenologically investigate the effect on primordial GWs spectrum induced by 
varying c^, we assume that the variation of c? can be described by a step-like function 


c T 


c T i (t < r 0 ) 

C-T2 0 > r 0 ) 


( 5 ) 


where To < 0 is the transition time. 

We take the background evolution to be the slow-roll inflation. Of course, it is also 
provided that the sudden change of ct won’t affect the background evolution. Then, we 
have z't!zt = a"/a ~ (2 + 3c)/t 2 , and the equation of motion (4) becomes 

u k + ( c rk 2 - u k = 0 , ( 6 ) 


where v = y/^ + 3e~| + e. The solution to Eq.(6) is familiar, we can write it as 

Uki = \/-crikr [C^H^i-CTikr) + C^ 2 H (2) (-CTikr)] , t < t 0 , 

u k2 = \Z-c T2 kT [C2,iH^\-c T 2kr) + C 2 , 2 H^ ) {-^kr)] , t > t 0 , (7) 


where and H^ are the first and second kind Hankel functions of z/-th order, respectively. 
These coefficients C are functions of the comoving wave number fc, but constants with respect 
to conformal time t. C± y i and Ci ,2 are determined by the initial condition. 

Here, we set the initial condition as the standard Bunch-Davies(BD) vacuum. Therefore, 
when crik 3> —, which corresponds to perturbations deep inside the horizon, 


Uk ~ 


_ e ~iCTikT 

\j2c T \k 


( 8 ) 
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When c.t\ k — , Uki in Eq.(7) should approximate to Eq.(8). Allow for the Hankel function 


v z 4 ) and Hy 2 \£) = J\e ^ I) when |£| —>■ 00 , we get 


r( 2 )rc\ — 


C 1A = 


IT 


C 12 — 0. 


( 9 ) 


2\/ orik 

The coefficients C 2j i and C 2i2 are determined by requiring and w/ to be continuous at 
r = tq, i.e. the matching condition. Then we obtain 


C r 2,i — 


C 22 — 


16VCT2 


16 y/CT2 


c T i yH^l +v (—c T ikr 0 ) - H^ u (-c T 1 kr 0 )j H® (-c T 2 kr 0 ) 

+ c T2 (-H {2 l +u (-c T2 kT 0 ) - H[ 2 ^(-CT 2 kT 0 )^ Hjp(-CTihT 0 ) 

T ° k C T l(^-H^l +l/ {-CTlkTo) + H^l^-CTlkTo^j H^\-CT2kTo) 

+ c T 2 (H-l +v (-CT2kT 0 ) - H^ui-^kro)^ H^\~c T ikTo) 


( 10 ) 


The spectrum of primordial GWs is defined by P T = (k 3 /2n 2 )(0\h i j(r, — k)hij(r, k)|0) 
with r —> 0, which is only a function of comoving wave number k. After neglecting the 
slow-roll parameter, from Eq.(7) with u = 3/2, we have 
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\ u k2\ ~ -7- 7 = |C*2,1 — C 2 , 2 | 

—C T2 kT V7T 


Therefore, we obtain the power spectrum of primordial GWs as 

Pt = 77 / Mt> fc)|2 = p Pq^4 2 |C ' 2 ’ 1 ■ C2 - 212 ■ 

where 


(ii) 


( 12 ) 


(13) 


PP = IHjnsK 

is that of standard slow-roll inflation without modified gravity, i.e. Qt = 1 and cti = ct 2 = 
1, and H in f is the Hubble parameter during inflation, which sets the scale of inflation. 

The effect of varying ct is encoded in C 2 j and (7 2)2 . We set x = ct 2 /cti and defined a 
new function 


f[k,k«,x)= i -^-\C v -C v 2 f, 
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where k 0 = —1 /(ctiXT 0 ). Then, the GWs spectrum (12) may be rewritten as 

-)inf f(k,ko,x) 


Pr = Pt 


Cj^iQt 


(14) 


(15) 
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FIG. 1: The function f(k,ko,x), where x = c^/cti- We have set x = 0.9 on the left panel and 
x = 1.1 on the right panel, respectively. 


where f{k, k 0 ,x) is obtained as 


f(k, ko, x) = -4 sin 2 (\ 


x 2 V k 0 J ' x 4 


( ^ \ /-i 2\ ^0 • / ^ \ 


(16) 


We plot f[k, k 0 ,x) with respect to k/k 0 in Fig.l, in which we set x = 0.9 on the left panel 
and x = 1.1 on the right panel, respectively. Here, the transition time To = —1 /(cti^A:o) 
sets a character scale l/fc 0 - When k ko, i.e. the perturbation mode has longer wavelength 
than 1/fco, we have f{k,ko,x) ~ 1, and Pt = P? /(Cj^Qt) is scale invariant, which is the 
result of slow-roll inflation with almost constant ct and Qt [23] [34], When k ko, we have 


f(k,ko,x) 


x A 


00 ^( 4 ) 


(17) 


thus f[k,ko,x) oscillates between 1/x 2 and l/x 4 , and Pt oscillates correspondingly, just as 
we can see from Fig.l. 

In the case of varying sound speed cs of primordial scalar perturbations, the sudden 
change of cs may lead to the oscillating modulation to the primordial scalar spectrum, as 
well as the CMB TT-mode spectrum, just as found in [29]. In addition, the oscillation 
in the primordial scalar spectrum can also be attributed to some other effects, such as 
inflaton potential with a small oscillation [35] [36], a sudden change in inflaton potential or 
its derivative, e.g., [37] [38] [39] [40]. Thus the oscillation in the primordial scalar spectrum 
may be implemented without modified gravity, as has been mentioned. 

However, the oscillation in the primordial GWs spectrum can only be attributed to the 
modified gravity. When the gravity is GR, Pt equals to P ^ 1 ^, which is given in Eq.(13). The 
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FIG. 2: Theoretical CMB BB and TT-mode power spectra for our oscillating GWs spectrum 
(15) (brown line in the left panel and red solid line in the right panel) and the power-law GWs 
spectrum for reference(blue dashed line in the left panel and green dashed line in the right panel). 
The insets of the right panels are the TT-mode spectra for our oscillating GWs spectrum (the 
yellow solid lines) and the power-law GWs spectrum (the blue dashed lines) for reference. We set 
r = 0.05 and ko = 1/30000 Mpc~b 

oscillation of P™‘ certainly requires H in f is oscillating, which is impossible, unless the null 
energy condition is violated periodically. Though the particle production may also modify 
the GWs spectrum [41] [42] [43], it only leads to a bump-like contribution, which is entirely 
different from the behavior of oscillation. 

III. CMB B-MODE POLARIZATION SPECTRUM 

The primordial GWs is imprint in CMB as the B-mode polarization. Thus the oscillation 
in the primordial GWs spectrum will inevitably affect the B-mode polarization spectrum. 
To see such effects, we plot the CMB BB and TT-mode correlations in Fig.2, in which 
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FIG. 3: BB-mode spectra at low multipoles for our oscillating GWs spectrum (15) with different 
x(solid lines) and the power-law GWs spectrum with different r r j(dashed lines). 

Pf’^ in (15) is parameterized as 

inf 

/l.\ n T 

P'r’ = rK 1 (j-j • (18) 

Here, we assume that the scalar spectrum is hardly affected by the modified gravity, which 
will be clarified in Appendix A. Thus the scalar perturbation spectrum is set as = 
A% f (k/k*) n K f ~ 1 , in which A’^ is the amplitude of the scalar perturbations. In addition, 
we also assume that after inflation the propagating speed ct is unity, so that the spectrum 
of B-mode polarization is not affected by relevant evolution at late time, or see [11],[12], 

In the left panels of Fig.2, we see some obvious enhancements or suppressions to the 
reionization bump in the BB-mode spectrum, which depend on the oscillating effect on 
corresponding scales. The height of the reionization bump can be estimated roughly as [44], 

» ^ 5 (l-e- r ") 2 C^ 2 , (19) 

where Cff stands for the TT-mode spectrum from the primordial GWs without the reioniza¬ 
tion and T r i is the optical depth to the beginning of reionization. The periodic enhancements 
and suppressions of the reionization bump are a reflection of the oscillations of primordial 
GWs spectrum on large scales. In addition, we can also see some obvious oscillations around 
the recombination peak at l ~ 80. 
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In the right panels of Fig.2, we see that the TT-mode spectrum is hardly affected by 
the oscillating primordial GWs power spectrum, since the contribution of GWs to TT- 
mode spectrum is negligible, compared with the scalar perturbations. The case of EE¬ 
mode polarization spectrum is actually also similar. Thus the main effect of varying speed 
of primordial GWs is on the B-mode polarization, which makes the B-rnode polarization 
spectrum show its obvious enhancements or suppressions to the reionization bump and 
oscillations around the recombination peak. 

In Ref. [45], the authors have pointed out that it is possible to set to one the propagating 
speed of GWs by a proper redefinition of the metric. They got the gravitational waves 
spectrum in their Eq.(12) “same” as the standard one in the form. However, since they have 
made a redefinition of the time coordinate and the scale factor in Eq.(9), the variation of 
their H with respect to t after redefinition is not the same as the variation of H with respect 
to t. Therefore, the oscillating feature induced by a step-like ct is encoded in H, the result 
in both frame should be the same. 

It might also be a concern whether such a B-mode polarization spectrum can be imitated 
by adjusting other cosmological parameters or not. In Eq.(19), the optical depth th is 
relevant with the height of the reionization bump. We show the BB-mode spectrum with 
different r r j in Fig. 3. We see that the change of r r ; can only alter the overall amplitude 
of BB-mode spectrum at low multipoles, but hardly create the oscillation. In addition, in 
inflationary models with pre-inflation era, the reionization bump could also be suppressed 
(e.g. [46] , [47] , [48] , [49] ) or enhanced (e.g. [50]). However, similarly, these models also only 
alter the overall amplitude of the BB-mode spectrum, without oscillation, at low multipoles. 
These results indicate that although the BB-mode spectrum may be modified by other ways, 
the oscillating modulation leaded by varying the speed of primordial GWs is difficult to be 
imitated. Thus the measure of B-mode polarization spectrum provides an appropriate way 
for testing the corresponding gravity physics in the primordial universe. 


IV. DISCUSSION 

In low-energy effective string theory and modified gravity theories, the propagating speed 
Ct of primordial GWs may deviate from unity. We calculated the spectrum of primordial 
GWs, assuming that ct has a step-like variation and the background of slow-roll inflation 
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is not affected by it. We found the spectrum of primordial GWs acquires an oscillating 
modulation, which makes the B-mode polarization spectrum show its obvious enhancements 
or suppressions to the reionization bump and oscillations around the recombination peak. 
The intensity of the modulation is determined by oril&ri- The frequency of the modulation 
is determined by k 0 = —1 /(ct2Tq)- Both depend on the dynamics of theoretical models. 

The oscillating behavior of the B-mode polarization can only be attributed to the effect of 
modified gravity, since it can hardly be imitated by adjusting other cosmological parameters. 
The oscillating behavior is so rich in feature that it may be easily discriminated from the 
variation of other parameters or other features, thus the upcoming CMB experiments, such as 
CMBPol, B-Pol, will be able to put a tight constraint on the propagating speed of primordial 
GWs, and so the corresponding theories, provided that the foreground contribution is under 
control. In a certain sense, our paper again highlights the significance of B-mode polarization 
measures in exploring the fundamental physics of primordial universe. 

Here, we only postulate a simple step-like variation of ct, which, however, might be far 
complicated in some modified gravity models, as well as accompanied by the variation of Qt- 
The effect could be non-trivial in more general cases, which is under study. But in a certain 
sense, a smooth change of ct will induce oscillations too, see e.g. [51] for the case of scalar 
perturbations, so we expect that the case of GWs is similar. When we focus on the B-mode 
polarization spectrum, the assumption we adopted is that after inflation the propagating 
speed ct is unity, which may also be relaxed. Moreover, it may well be possible to produce 
oscillatory features beyond the standard slow-roll background. The varying ct and Qt will 
also affect the non-Gaussianities of primordial perturbations. The relevant issues are open. 
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Appendix A: 

In this Appendix, we will argue how to realize the change of ct and the changeless of Qt 
in (1) in low-energy effective string theory and modified gravity theories. 

In low-energy effective action of string theory, the simplest extension of the lowest-order 
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action is e.g.[17] 


£correction ~- ^ ~ ( C 1 ^GB + C 2 G*™ 3 pSpdvty) ( 20 ) 

where G tw = R^ v — g^R/2 , and R 2 GB = R. l ,, v \ P R! ll ' Xp — 4 R^R^ + R 2 is the Gauss-Bonnet 
term, a' is the inverse string tension, A is a parameter allowing for different species of string 
theories, c\ and C 2 are coefficients. We have neglected the terms with Utp and (d^d^tp) 2 , 
since both do not contribute to GWs. 

The introducing of (20) will affect not only the GWs, but also the adiabatic scalar pertur¬ 
bations, of course, it is interesting to check its effect on the latter. However, for our purpose, 
we will regard as the spectator field, so that the effects of (20) on the background and the 
scalar perturbations are negligible. 

The action for GWs is (1), and [17] 


Qt 


c 


2 

T 



(%ciiH inf - C 2 £0 2 ) , 
(c 2 £0 2 + 8c, ^ 


( 21 ) 


see also Ref. [52] for that with c\ = 0. When 8ci£H in f = c 2 £yj 2 is imposed, Qt = 1. Then, 
we have 

Ct = 1 — 4cia'A(£Hi n f + £). ( 22 ) 


The step-like variation of requires that £,H in f + £ has the step at r = r 0 . 

As an example, we will give a numerical result of the variation of c\ and <p with respect 
to cosmological time t in the case of x — ct 2 /cti = 10. According to Refs. [14] [17], we adopt 


£((p) = —e Ci = —1, a' = 1, and A 


1 

8 


(for Heterotic string). 


(23) 


We take the expression of ip(t) as 
ip(t) — t — In 


Ae t ~ t 1 

ffie 4 + b 2 A — He*- -j= —I- Ate 1 + Ae 2 erf(- t) — Ae t (t — l)erf(t) 

7r 2 


. (24) 


where A = —1 + l/x, and erf(i) = f* e~ f2 dt is the Gauss error function, b\ and b 2 are 
some constants. We plot ip and c\ with respect to cosmological time t in Fig.4. 

Of course, the variation of could be more complicated than a step-like evolution, which 
would be harder to deal with. However, as it may, Qt > 0 and > 0 should be required 
to avoid ghost and gradient instabilities. 


11 









FIG. 4: The variation of ip and c|, with respect to cosmological time t in the case of x = 0 ^ 2 / c ti = 
10, i.e., A = —0.9, where t = 0 corresponds to to- We have set b\ = 50, 62 = —50. 

The action (20) is actually equivalent to a subclass of the Horndeski theory, see Ref. [23], 
and so the analysis is also similar for the Horndeski theory. 
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